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The HIRFL-CSR external-target experiment (CEE) under construction is expected to provide 
novel opportunities to the studies of the thermodynamic properties, namely the equation of state of 
nuclear matter (nEOS) with heavy ion collisions at a few hundreds MeV/u beam energies. Based 
on Geant 4 packages, the fast simulations of the detector responses to the collision events generated 
using transport model are conducted. The overall performance of CEE, including spatial resolution 
of hits, momentum resolution of tracks and particle identification ability has been investigated. 
Various observables proposed to probe the nEOS, such as the production of light clusters, t/ 3He 
yield ratio, the radial flow, m~ /r* yield ratio and the neutral kaon yields, have been reconstructed. 
The feasibility of studying nEOS beyond the saturation density via the aforementioned observables 


to be measured with CEE has been demonstrated. 


I. INTRODUCTION 


The phase diagram of nuclear matter is of significance 
in understanding the properties of the strong interaction 
described by Quantum Chromodynamics (QCD) and the 
evolution of the universe. With the relativistic heavy ion 
experiments, the various experimental signals have re- 
vealed that strongly interacting matter may undergo the 
transition from the hadron phase to the quark-gluon- 
plasma (QGP) phase, for which the transition at low 
baryon density and high temperature is crossover [1-3]. 
Fig.1 shows the QCD phase diagram of current under- 
standing [4]. Despite of great progress along this di- 
rection, a few important questions remain unsolved con- 
cerning the QCD phase diagram. For instance, what are 
the thermodynamic properties of QCD and where is the 
phase boundary between QGP and hadronic matter at 
high net-baryon density? Does the critical end point 
(CEP) exist and where is the location? To shed light 
on these questions, it is of great significance to explore 
further the QCD phase diagram with more systematic 
studies. 

Relativistic heavy ion collisions (HIC) over wide en- 
ergy range provide the unique way to study QCD phase 
diagram in terrestrial laboratory. Over the past two 
decades, nuclear collisions at the Relativistic Heavy Ion 
Collider (RHIC) and the Large Hadron Collider (LHC) 
have collected a large amount of data searching for first- 52 
order phase transition signals and precise phase transi- 53 
tion CEP. Between 2010 and 2017, RHIC completed the 54 
first phase of the beam energy scanning program (BES- ss 
I) [5, 6]. The second phase (BES-II) was launched in se 
2018, focusing on the collisions with VS < 27 GeV tos 
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FIG. 1. (Color online) QCD phase structure diagram of 


temperature and baryon chemical potential description [4]. 
The red line in the figure is the chemical freeze-out curve of 
common nuclear substances obtained by experiments. The 
black square points and black solid lines represent the ex- 
pected first-order phase boundaries and critical points. The 
red circle represents the CEP at which the nuclear material 
becomes gas-liquid phase. The study areas of RHIC beam en- 
ergy scanning and fixed target (USA), NICA (Russia), FAIR 
(Germany), and HIAF (China) are shown at the top of the 
figure. 


search sensitive observable of QCD phase transitions [7- 
9], including net proton [4, 10-13], net charge [14], net 
Kaon multiplicity distribution and high order moment 
analysis [15]. Fig.2 shows the energy dependence of net 
proton high momentum «ko? in Au+Au collisions with 
the centrality of 0 — 5% collected in RHIC BES-I [4]. It 
is observed that ko? of net proton is closely related to 
that of proton and exhibits a non-monotonic dependence 
on collision energy with a significance of 3.1o [11, 13]. 
In addition, the prediction of UrQMD calculation, which 
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incorporates conservation laws and most of the relevant ss 
physics apart from a phase transition and mean field po- s7 
tentials, does not coincide with the trend of the data ss 
points at low energies, which means that accurate mea- ,, 
surement is very important in the supra-saturation den- ,, 
sity region and lower beam energies to further reduce ex- ,, 
perimental errors. In addition to the existing ones, many ,, 
more experiments are in schedule, like CBM at FAIR, ,, 
MPD at NICA etc. zi 
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FIG. 2. (Color online) Energy dependence of net proton and™? 
proton Ko? in Au+Au collisions at the center (0—5%) of RHIC! 
BES-1 [4]. The results of the HRG and UrQMD models are 
shown as dotted lines and gold bands, respectively. 115 
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In hadron phase, the key component of the QCD phase” 
structure is the equation of state of nucleonic matter”? 
(nEOS) [16]. For the convenient definition of nEOS,™ 
one writes the nucleon specific energy E(p,6) of nucle- 


onic matter at zero temperature in quadratic form a 
122 


123 
( 1 ) 124 


where 6 = (fn — Pp)/(Pn + Pp) is the isospin asymme- 
try, Pn and pp are the density of neutrons and protons. 
Eo(p) is the nucleon specific energy in symmetric nuclear 
matter and Esym(p) is the density dependent nuclear 
symmetry energy. Keeping the lowest order terms, for 
the symmetric nuclear matter, the nEOS can be written 


in 

Jo Ca _ 20 (= a 

6 \ 30 "24 \ 390 : 

(2) 

where the coefficients xo, Jo and Zo are the incompress- 
ibility, skewness and kurtosis parameters, respectively. 


The nuclear symmetry energy is expressed as 


E(p, ô) = Eo(p) Eg Exym(p)6" 


2 
as Pm 
Eo(p) = Eo(po)+ o (e e) i 

2 3p0 


Esym(p) = Esym(po) 
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where L and Ksym are the slope and curvature param- 
eters of the nuclear symmetry energy at normal nuclear 
matter density po, respectively. 


In HIC from a few hundred MeV/u to a few GeV/u, 
nuclear matter above saturation density is created. The 
density gradients of the compressed nuclear matter will 
be developed to various forms of flow carrying the infor- 
mation of nEOS, and in return, the flow has been used to 
probe the incompressibility of nuclear matter [17]. Re- 
cently, HADES collaboration has published the proton 
triangular flow v3 data at 2.4 GeV. Compared to UrQMD 
transport model calculations, it is shown that the proton 
vg exhibits high sensitivity to the hadron medium EoS 
[18]. In addition to the flow analysis, the production of 
K meson, particularly near threshold, is sensitive to the 
density of the compressed participant region, and thus 
is employed as a probe of nEOS. Thanks to the great 
efforts from the communities of nuclear physics and as- 
trophysics, the knowledge of Eo (p) from po to a few times 
of po is accumulated. 


The least known term of nEOS is the density behavior 
of symmetry energy Esym(p). Since the Esym(p) gives 
rise to the isospin diffusion and isospin drift mecha- 
nisms, in the terrestrial laboratory, HIC involving neu- 
tron rich nuclei provide unique opportunities to constrain 
Esym(p) near and above the saturation density. In a wide 
energy range, various isospin probes have been identi- 
fied and applied to constrain Esym(p) near and above 
po in the last two decades [19-26]. Some new probes 
are also introduced or updated very recently, including 
the neutron skin thickness of heavy nuclei [27, 28], the 
angular distribution of average N/Z of light charged par- 
ticles [29], anticorrelation of the t/*He yield ratio with 
the heavy clusters [30] for p ~ po and the 7~/zt ratio 
for p > po [81]. Moreover, since the observation of the 
neutron star merging event GW170817, great progress 
has been made to constrain Esym(p) by combining the 
GW170817 and heavy ion data [32, 33]. For a recent 
review, one can refer to [34]. 
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FIG. 3. (Color online) Research the constraints of symmetry 
energy in the supra-saturation density region from various 
theories, and display the constraints on the boundary through 
neutron star data [35] 
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Despite of great progress achieved in the last tworso 
decades, still the constraint of Esym(p) suffers from large 
uncertainties [33, 36-39], particularly at high density,,, 
region. World wide efforts in this direction is still on- 
going. In order to finally determine Ezym(p) in the, 
supra-saturation density region, some new experiments 
are scheduled or have conducted at some large scientific y 
devices in the world, including but not limited to HADES 
at GSI, StRIT at RIKEN, RAON-LAMPS in Korea etc. s 
These latest experiments will certainly provide better eX- a 
perimental conditions and improve the statistics and Te- as 
duce the uncertainties. For the recent review, one Can yo 
refer to the long report [40]. 
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191 


To this end, in addition to the RHIC-STAR energy” 
scan experiment, there are several large experimental fa-’” 
cilities in operation or construction around the world,’ 
which are of great significance for mapping the QCD” 
phase diagrams at high net baryon number density, as 
well as for constraining the nEOS in hadron phase. In ad-"” 
dition to the aforementioned facilities around the world,” 
CEE is currently under construction at the T4 Exper-” 
imental Terminal of the Heavy Ion Research Facility” 
(HIRFL) in Lanzhou [41, 42], aiming at collecting the™ 
first physics data by 2025. The High Intensity Heavy” 
Ion Accelerator Facility (HIAF) is under construction in” 
Huizhou, China [43]. The HIAF in design will be able™ 
to provide beams of whole ion species, with the highest?” 
center of mass energy in U+U collisions being V SNyN = 4 
GeV. The upgraded CEE?* will be placed in terminal T2 
of the HIAF for high energy experiments, providing great 
opportunities in the studies of QCD phase structure in 
high baryon density region and EoS of nuclear matter in 
high density. 


In this paper, the design, performance and feasibility 
of various physical goals related to the EOS studies with 
CEE are introduced. The simulations are carried out in 
the framework of Geant 4. The ultra-relativistic quan- 
tum molecular dynamics (UrQMD) model which incorpo- 
rates the Skyrme potential energy density functional and 
the in-medium nucleon-nucleon cross section is adopted 
as the event generator [44-46]. The analysis and re- 
constructions of the simulated data are done in CEE- 
ROOT platform, which is developed by the collaboration 
of the CEE experiment, inheriting from Fair ROOT [47]. 
The paper is arranged as following. Section II presents 
the structure of CEE in the technique design and the 
main functions of each detector subsystem. Section III 
introduce the expected performance of CEE, including 
the spatial acceptance, momentum resolution, PID abil- 
ity, and the construction of CEE global trigger signal. 
Section IV presents the verification of a few probes sen- 
sitive to the EOS of nuclear matter beyond saturation 
density, including the yield of light clusters, radial flow, 
a /x* ratio and yield of neutral Kaon mesons. Section 
V presents a summary. 206 


II. STRUCTURE AND SIMULATION OF CEE 
A. Structure of CEE 


Aiming at the studies of the EOS of nuclear matter 
at approximately 29, where po is the saturation den- 
sity, CEE is designed to measure the charged particles 
in nearly 47 space in center of mass reference [41]. The 
design of the CEE is shown in Fig.4. The main compo- 
nent of CEE is a large-gap dipole magnet, housing the 
tracking detectors in a nearly homogeneous field of 0.5 T 
along the vertical direction. The target is located inside 
the field, so that the tracking detectors can cover large 
solid angle. The tracking detectors include two time pro- 
jection chambers (TPC) [48] surrounded by the inner 
time-of-flight (ITOF) detectors covering the midrapid- 
ity [49], and an array consisting of three multiwire drift 
chambers (MWDC) followed by an end cap time-of-flight 
(eTOF) wall for the measurement of charged particles at 
forward rapidity [50]. At the downstream end of the 
CEE spectrometer, there is an array of plastic scintilla- 
tors to measure the charged particles near zero degree 
(ZDC) [51]. The start timing detector (To) [52, 53] and 
the active collimator (AC) are placed on the beam line in 
the upstream side of the target. Besides, there is also a 
silicon pixel positioning detector (SiPix) [54, 55] to mon- 
itor the beam position at the upstream side to Tp. Tab.I 
shows the main technical performance of CEE system. 


@ TPC 


© iTOF 
© Dipole magnet MWDC 
© AC & SiPix © eTOF 
© To zDC 


FIG. 4. (Color online) The design of the CEE 


TABLE I. (Color online) Technical indicators of CEE. 


Item value 
Maximum beam energy 0.5GeV/u(U) — 2.8GeV(p) 
Bean type prt 
Maximum event rate 10 kHz 
Acceptance > 50% 
Total channel number 20k 


The tracking detector at midrapidity are the two iden- 
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tical TPCs surrounded by iTOF walls. The two TPCs arezes 
installed in a left-right symmetric mode in order to leavezss 
the 20 cm space in between for the beam passing through.2.7 
The target position is in the gap of the two TPCs, beingzss 
15 cm to the rear surface of the field cage. The electronicsz«s 
of the TPCs are placed on the top plane above the two27 
TPCs. The high voltage is fed into the field cage frome 
the bottom plane. A high voltage step-down structurezz. 
surrounded by teflon blocks is designed to avoid a sharpz73 
voltage decrease near the HV point. The field cages arez 
made of Kapton layer with copper electrodes printed ons 
both surfaces, ensuring the distortion of the electric field. 
to be less than 0.1%. A three-layer large-area gaseouso77 
electron multiplier (GEM) foils are stacked on the top of» 
the sensitive volume of TPC to amplify the arrived elec-27 
trons and then induce signals on the readout pads. Thezso 
signals of each pad are readout by the specially designedzsı 
front end electronics (FEE) and then transferred to ther 
digitization board adopting SAMPA chips [56, 57]. 283 


The forward tracking detectors consist of three? 
MWDGs and an eTOF wall placed at forward angle. Two? 
sets of MWDC are placed inside the field of the dipole 
and the third one followed by the eTOF wall are placed? 
outside the magnet. On each MWDC and the eTOF wall, 
an inactive area is designed for the beam particles to pass 
through without inducing enormous charges in the sen- 
sitive volume of the detectors. Each MWDC contains 
six layers of sensitive wires in three directions, called X, 
U and V wires, meeting 0°, 30° and —30° with respect 
to the vertical direction, respectively. For each direction 
of wire, there are two layers of drift cells, displacing by 
half size of a cell in order to discriminate the ghost hit. 
Both the entrance and exit windows of the MWDCs are 
made of aluminum coated Mylar foil. The detector op- 
erates in atmospheric flow-gas mode. The high voltage 
are fed to the anode wires directly. The signals, readout 
with a capacitance coupling, is amplified and shaped in 
the FEE before being transferred to a SCA chip for dig- 
itization. Then the digitized data are sent to the BDM 
module, where the timing and amplitude information is 
online computed and transferred to DAQ system for saV- 
ing. Tab.II lists the main geometric dimensions and tech-,s 
nical specifications of the TPC and MWDC. 566 


The time-of-flight (TOF) signals are measured by To, 
inner TOF wall (iTOF) and the outer TOF wall (eTOF),22 
respectively. The To placed in front of the target is29 
made of scintillator foil read out by SiPM from the side,294 
delivering the initial time of heavy ion collision. Bothzss 
iTOF and eTOF, consisting of multi-gap resistive platezs 
chambers (MRPC), are installed to record the arrivals 
time of the charged particles passing through TPC andzs 
MWDC, respectively. Thus, combining the tracking and2% 
TOF detectors, one can deduce the PID information viasoo 
dE/dx — p/Q in the tracking detectors and p/Q — 6,30 
respectively, where p is the momentum and £ is the ve-s02 
locity calculated from the track length and TOF. Tab. IIIs 
shows the main performances of TOF system. The totalso 
area of the iTOF and eTOF systems is about 12 m?. Inso 


iTOF and eTOF, the hitting position can also be derived 
through the time difference of the signals at both ends. 

The iTOF, grouped by modules, surrounds the left, 
right and bottom surfaces of the TPC. The upper sur- 
face is not covered by TOF detectors because the read 
out electronics of the TPC is arranged on the top. The 
left and right iTOF walls are composed of six iTOF mod- 
ules, each iTOF module contains three MRPC detectors, 
and the left and right envelope size is 100(x) x 1000(y) 
x 1355(z) mm’. The bottom iTOF wall contains two 
symmetrical parts, consisting of 6 iTOF modules and 12 
MRPC detectors, with the bottom envelope size of 670(x) 
x 1355(z) x 100(y) mm3. The eTOF wall is located at 
the down stream end of the dipole magnet, being 2.6 m to 
the target position. Each MRPC has 10 x 0.25mm air 
gaps and self-sealing technology is applied for the first 
time. The whole wall is divided into 7 modules, each 
module contains 3 ~ 4 self-sealed MRPC [50]. The sig- 
nals of the TOF detectors are collected and amplified by 
the MRPC special electronics chip NINO [58] in FEE, 
and then the high precision time information is obtained 
through the high performance time digitization module 
(TDM) [59]. 
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FIG. 5. (Color online) The magnetic field distribution in the 
center of the secondary magnet. The left (a) shows the mag- 
netic field distribution in the X-Z plane, and the magnetic 
field shield tube is shown in the red box; The right figure (b) 
shows the distribution of magnetic field intensity in the y di- 
rection in the beam direction 


The active collimator functioning to veto the off-target 
reaction background is placed between To and the tar- 
get. It is composed by four plastic scintillators ( 250 x 
250 mm?) distributed around the entrance of the target 
chamber, leaving a 50 x 50 mm? hole in the center. The 
SiPM is used to read out the signals. 

The ZDC is mainly used to measure the spectators 
from HIC, and to determine the centrality and the reac- 
tion plane of collisions. The ZDC is located at the most 
downstream position of the whole spectrometer, about 
300 cm to the target. The overall transverse layout of 
ZDC is rotationally symmetric to the beam line. It is 
divided into 24 equal sectors, each of which contains 8 
trapezoid scintillators placed in row along the radial di- 
rection. The outer radius of the ZDC sensitive area is 
100 cm while the inner radius is about 10 cm, through 
which the non-reacted beam particles can pass. 

The large gap superconducting dipole magnet is to 
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TABLE II. (Color online) The geometry size and technical indicators of TPC and MWDC 


TPC parameters 


MWDC parameters 


Item value 


Item value 


Effective sensitive area (half) 


500(x) x 800(y) x 900(z) mm3 


MWDC1 effective sensitive area 136(x) x 63(y) cm? 


Channel number > 10000 MWDC2 effective sensitive area 190(x) x 88(y) cm? 
Sensitive area of GEM film 900 mm x 500 mm MWDCS effective sensitive area 274(x) x 126(y) cm? 
Track multiplicity < 100 Channel number > 3000 
Maximum event rate 1kHz Transverse position resolution 0.3 mm 
Maximum field strength 400 V/cm Track reconstruction efficiency ~ 98% 
Momentum resolution ~ 5% Momentum resolution ~ 5% 
TABLE III. (Color online) The technical indicators of iTOF and eTOF 
iTOF parameters eTOF parameters 
Item value | Item value 
Total effective area 6 m? Total effective area 3.2 x 1.6 m? 
Channel number ~ 2400 Channel number ~ 1500 
Average count rate of single channel 10kHz Average count rate of single channel 10kHz 
Time resolution accuracy < 50ps Time resolution accuracy < 60ps 


bend the charged particles produced in the projectile-sss 
target collisions. The raw dimension of the magnet. isss 
2700(L) x 4300 (W) x 2.7 (H) mm8. The space of thes 
uniform field housing the sensitive volume of the TPC issa 
900 (L) x 1200 (W) ~ 800 (H) mm°. The central fields: 
is 0.5 T in strength with a uniformity better than 5%.s43 
Fig.5 shows the magnetic field distribution in the centerass 
of the dipole magnet. 345 

In addition, the CEE also includes the following sup-s4s 
porting systems. The data acquisition system (DAQ): 
uses D-Matrix architecture to provide data storage. Thess 
main function of the clock system is to provide high preci- 
sion global synchronous clock for the electronics modulesaso 
and DAQ system. The trigger system provides globalss: 
trigger signals in beam experiments or in detector testas2 
run. Refer to Section IIIC for details about the con-ss3 
struction of the global trigger signals of CEE. 354 
355 
356 
B. Event generator and simulation tool 357 
358 


In order to constrain nEOS, particularly for supra-*” 
saturation density, one usually use the given Eo(p) and*° 
Esym(p) as input in transport model for HIC or in solv-** 
ing the Tolman—Oppenheimer-Volkoff (TOV) equations” 
for neutron star structures. The predicted results are 
then compared to experimental data or astrophysical ob- 
servations. By adjusting the parameters of Eo(p) and 
Esym(p) until the results converges to the experimental 
results, one arrives a constraint of nEOS. And in return, 
transport model is usually employed as the event gener- 
ator in the feasibility studies of nEOS via in terrestrial 
heavy ion experiment. 

UrQMD model is a transport model widely used to 
describe the relativistic heavy ion collisions and proton 


363 


induced reaction process [60, 61], ranging from sub- 
GeV/u to TeV/u energies covered by SIS, AGS, SPS, 
RHIC and LHC. In order to apply the UrQMD model 
as an event generator for CEE, the model is adopted by 
improving the isospin-dependent nuclear potential and 
collision terms. 

The UrQMD model as the event generator is operated 
and provided by Huzhou University in China [44-46], 
and the main improvements include: 1) The mean-field 
potential derived from the Skyrme potential energy den- 
sity functional has been implemented into [46, 62, 63] 
2) The nuclear medium effects on NN cross sections 
are considered, the density-, momentum- and isospin- 
dependent NN cross section is adopted [44, 64-66]; 
3) The isospin-dependent minimum spanning tree (iso- 
MST) method [67] is adopted to construct clusters. In 
addition, the yield and flow observables about a and 
K, which are effective and sensitive probes of symme- 
try energy, can be better reproduced by introducing 
proper pion-nucleon potential and kaon-nucleon potential 
[33, 68, 69]. The effective two-body interaction potential 
energy U is decomposed into the Coulomb energy Ucou, 
the Skyrme potential energy U, and the momentum- 
dependent potential energy Uma. The U, and Uma is 
written as U, ma = f Up,madr, where 


2 n+l 
Up 2A P E P 7 Isur (Vp)? 
2p n+l po 2po (4) 
Jsur,iso 2 sym 
HEE Y (Pn Hum, 
a [V (Pn — Pp) + up 


and 


L (5) 


Umd = tma In? [1 + ama(P; = p;)’] p 
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The symmetry potential energy density functionalss 
reads as 


398 


P P P \5/3| 52 
Asym(— + Bsym(—)? + Ceym(— Ô D. 399 

y nr y we, y a P 
(6) 


For the Skyrme interactions, the symmetry potential 
energy terms come from two-body, three-body interac- s 
tion terms, and its related parameters are Asym, Bsym, 
and Csym, respectively. All the parameters used in the 
UrQMD model, such as a, 8, N, gsur, Jsur,iso; Asym, Bsym, 
Csym can be derived from the standard Skyrme param- 
eters, £o, Z1, T2, T3, to, ti, ta, t3, Oo (62, 63, 70]. In ad- 
dition, one can set Asym = Csym = 0 and Bsym = C,/2, 
and adopt the density power law form to investigate the1os 
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density-dependent symmetry energy|70, 71] 406 
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The simulations and the reconstructions of CEE ex-”” 


periment is implemented in the framework called CEE-*” 
ROOT, a platform inheriting from Fair ROOT software” 
[47] based on the ROOT and Geant 4 packages. User“ 
routines are developed to complete the data calibration,” 
event reconstruction, efficiency correction and histogram-*” 
ing. In CEEROOT, the subdetectors, the dipoles and*” 
the main materials are constructed using Geant 4 class‘ 
libraries. When the particles produced in HIC provided*” 
by the event generator are input to the simulator, the‘ 
propagation of the particles and the responses of the de-”* 
tectors are simulated. Then the signals in the sensitive’? 
volume of the detectors are digitized and saved to hard” 
disk for further reconstruction. In this studies, the de-*” 
tector systems, including the magnetic field, To, iTOF,”* 
eTOF, TPC, MWDC and ZDC are constructed as de-*” 
signed. The improved UrQMD model described in Sec-*” 
tion.II A is used as the event generator. Fig.6 shows the” 
event display and detector configuration in the simula-*” 
tions. The detector performance based on fast simula-*” 
tions will be introduced in the next section. a 
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FIG. 6. (Color online) The detector structure and the event 
display of CEE in fast simulation. 
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II. CEE PERFORMANCE STUDY 


This section presents the fast simulation studies 
on CEE performance. Using a symmetric collision 
208 Pþ+208Pþp, the acceptance, the PID ability and the 
centrality selectivity are discussed. The design of the 
global trigger design in heavy ion collision of CEE is in- 
troduced. 


A. Acceptance 


It is usually required to measure nearly the whole prod- 
ucts in the studies of nuclear EOS and QCD phase struc- 
ture, naturally one favors large acceptance for CEE. Fig.7 
shows the phase space distribution of different kinds of 
particles at b = 0 fm in ?°°Pb+?°8Pb collision at 0.3 
GeV/u. The abscissa is the reduced rapidity defined by 
yO = Yab/Yem — 1, where yap is the rapidity of par- 
ticles in the laboratory reference, and Yem is the rapid- 
ity of the center of mass of the collision system. The 
ordinate corresponds to the reduced transverse velocity 
u® = p,/MByyv, where p, and m represent the trans- 
verse momentum and the mass of the particles, respec- 
tively. Here 6, and y are the beam velocity and the 
corresponding Lorentz factor in laboratory. The curves 
corresponding to different angles in laboratory reference 
are plotted on top of the data. The MWDC array cov- 
ers the range below 30°, while the TPC mainly covers 
iab > 30°. It is clear that since the target is placed in 
the space between the two TPC with a distance of 15 
cm to the rear end of TPCs, the space at large angles 
with iab > 90° are also partly covered. From the phase 
space distribution of protons, it is also shown that in the 
vicinity of 30° both TPC and MWDC have the detection 
ability and hence this region is overlapped. It is shown 
that more than 90% of the phase space in center of mass 
can be covered, particularly for the symmetric systems 
because one can reflect the distribution with respect to 
y®) = 0 following the physical requirement. 

Fig.8 shows the scattering plot of the momenta p vs. 
Mab for protons, tritium, t+ and m7 in the reaction of 
208PH+298Pb at 0.3 GeV/u in laboratory reference sys- 
tem. The distribution only pass the filter of TPC and 
MWDC with a track length cut Lirk > 30 cm, while 
no conditions of the TOF hits are applied. Again, the 
wide angular acceptance of the tracking detectors are il- 
lustrated. 


Fig.9 further presents the angular coverage of the 
tracking detector TPC and MWDC on 0 — ¢ plane for 
protons as an example. Particles with a track length 
shorter than 30 cm in the TPC have been cut. While the 
polar angle is well covered at iab < 100°, there is obvi- 
ous efficiency lose in azimuth. Namely in the vicinity of 
o ~ 90° and ġ ~ 270°, both detectors lose the detection 
ability. 
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FIG. 7. (Color online) The phase space distribution of variousse2 
produced particles in the 27°°Pb+?°%Pb collision at 0.3 GeV /t.s6s 
464 


465 


466 
4 02 467 
468 
469 
10 470 
471 
472 


1 473 


474 


475 


101 


ar reile rat 
100 140 180 


FIG. 8. (Color online) Correlation between the momentum 
p and polar angle ĝia» in laboratory for protons, tritium, m” 
and 7 in the reactions of ?°°Pb+?°8Pb at 0.3 GeV/u. 


B. Particle identification 


Since CEE has two tracking detectors, TPC andasse 
MWDC, and TOF detectors, both methods usings7 
dE/dx —p/Q and p— TOF can be applied, respectively..s 
In the simulation studies, for the dE/dx — p/Q correla-a7 
tion obtained by TPC and MWDC, the total energy lossaso 
dE is accumulated in each step when a particle is prop-ss1 
agating in the sensitive volume of TPC or MWDC filledas2 
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FIG. 9. (Color online) Angular coverage of proton by the 
tracking detectors TPC (a) and MWDC (b). 


with the working gas of (90% Ar + 10% CO2). When a 
complete cell is processed, the energy loss rate of dE/dx 
is calculated by summing the energy losses in all steps 
and smeared by a given variance. Since the distribution 
of dE/dx follows Landau formula showing divergence at 
large value, usually a high energy cut is applied to fill the 
histogram. The momentum p is obtained by fitting the 
hits generated by the ionization of the energetic particles 
in propagation. Tracking finding is not yet implemented, 
so the fitting is done directly to each individual track and 
a smearing is introduced directly according to the posi- 
tion resolution of the tracking detector. Fig. 10 presents 
the dE /dx—p/Q for TPC and MWDC, respectively. It is 
shown that different isotopes of Z = 1 and Z = 2 can be 
well identified in TPC. While in MWDC, due to the total 
sampling times is only 18, corresponding to 18 layers of 
anode wire planes, the determination of the momentum 
is moderately degraded, although the tracking residue is 
controlled in 300 um. 
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FIG. 10. (Color online) The predictions of particle identifi- 
cation in TPC and MWDC using the energy loss dE/dx as a 
function of magnetic rigidity p/q. 


With increasing of p/q, the separation between neigh- 
bouring bands on the dE/dx — p/q is reduced, particu- 
larly between pions and protons beyond 1 GeV/c. Thus, 
p—TOF correlation is applied to improve the PID ability. 
Once the track length Lirk is determined by the tracking 
detectors, the speed of the particle can be calculated and 
the mass m can be resolved at a given momentum p by 
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c2t2 
gi c L? : 
where c is the speed of light. This method can 
be applied for both configurations of TPC+iTOF and 
MWDC+eTOF, respectively. To achieve the intended 
performance, the timing resolutions of iTOF and eTOF 
are designed as 50 ps and 60 ps, respectively. Fig.11 
presents the velocity of the particles as a function of the 
momentum p measured in TPC. It is shown that the pi-” 
ons and Z = 1 isotopes are well identified. <n 
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FIG. 11. (Color online) The change of charged particle flight 
velocity with momentum measured by iTOF 


Fig.12(a) shows the mass and momentum correlation 
by combining the eTOF delivering the time of flight and 
MWDC providing the momentum information. The mass 
distribution is further plotted in Fig.12(b). It can be seen 
that at the beam energies of interest, the light charged 
particles with Z = 1 can be separated clearly, although a 
small distortion is occurred at the low momentum border 
due to the influence of the gap of the detectors, which is 
not corrected in the track length estimation. 
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FIG. 12. (Color online) (a) Mass m and momentum p scatter-s28 
ing plot obtained by eTOF and MWDC, (b) mass distributionsz 
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The resolution of the momentum reconstruction of thess2 


8 


tracks in B field is determined by the hit position resolu- 
tion and the track length L in the bending area, written 
as 


ee a ee 0.0523 | ; 
= + (cot 6 9 
D =N 03BL2,) Iam T Ctre)" (9) 


where B is the magnetic field strength and @ is the 
polar angle. N is the number of samplings, ox, and og 
are the hit position resolution in x-z plane perpendicular 
to the field B and the angular resolution, respectively. 
Xo = 1.2 m is the radiation length. Given the hit posi- 
tion uncertainty of the tracks in the tracking detectors, 
the resolution of the momentum reconstruction can be 
investigated by fast simulations. 

Fig.13 shows the resolution of the momentum recon- 
structed from the hit information collected from sensitive 
volume of TPC for 0.3 GeV /u 2°8Pb+?°8Pb at b = 0 fm. 
Oxz = 0.6 mm and og = 0.5° are applied. It is shown 
that the typical resolution of 7+, proton and deuterium 
is about 5%. With the increase of momentum or in the 
forward angle region, the momentum resolution is grad- 
ually degraded, suggesting the necessity of using track- 
ing detector array at forward angle for the dipole-type 
spectrometer. Namely, the MWDC array is designed for 
CEE. The degrades of the resolution at low momenta is 
mainly due to the loss of tracking efficiency because the 
tracks with low momentum leave very short length in the 
sensitive volume. 
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FIG. 13. (Color online) Momentum resolution of TPC in fast 
simulation 


C. 


Global trigger signal construction 


In the beam experiment of the CEE, the construction 
of the global physical trigger signal has a twofold goals: 
1) First, one has to select the collisions on the target and 
suppress the background events on the beam path and 
2) It can provide trigger signals for the events at given 
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geometry, usually called as collision centrality. Similarse 
scheme is usually adapted in the existing spectrometers, 
for instance, as by FOPI [73] and SAMURAI [74]. 
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FIG. 14. (Color online) 0.5 GeV/u ?°8Pb+°8Pb reactions 
system, under the collision parameter b=0,1,2,3,4,5,6fm, these+ 
multiplicity of charged particles changes with the collision pa-sss 
rameters before and after passing through the detector. Thesss 
red (black) curve is the data before (after) the detector filter-sg7 
ing. pT 

At CEE, the multiplicity of charged particles firing the 
fast detectors, including iTOF and eTOF, is employed to 
construct the global trigger because it can provide an ap- 
proximate measure of the collision centrality. Figure 14 
presents the multiplicity of charged particles as a function 
of the impact parameter b in the transport model simu- 
lations of the reactions ?°°Pb+?°8Pb at 0.5 GeV/u. The 
multiplicity distributions before (open) and after (solid) 
the filter of the iTOF and eTOF are presented. It can be 
clearly seen that the multiplicity of charged particles fir- 
ing the iTOF and eTOF is anti correlated with b, giving 
certain ability to select the event centrality as the trigger 
signal of the CEE spectrometer. 


IV. STUDY ON THE EOS OF NUCLEAR 


MATTER 


In this section, we present the feasibility of studying 
the nuclear equation of state using several selected ob- 
servables at CEE. In the HIC in the energy domain of 
hundreds of MeV/u, the nuclear matter can be com- 
pressed to approximately 299, where pọ is the saturationsss 
density. In addition, due to the large stopping and largesso 
spatial-temporal volume, the sensitivity of the observablese1 
on nEOS is expected to be enhanced [75, 76]. Thus, itso 
is of interest to check whether the information of EOSse3 
carried by various the observables is discernible at CEE.so4 


A. light nuclear production 


One of the physical goal of CEE is to measure sys- 
tematically the production of 7H and *He at the beam 
energies between 300 and 600 MeV/u. The motivation 
of this measurement is to understand the correlation be- 
tween clustering and the isospin transport in HIC before 
arriving at a more stringent constraint on nuclear sym- 
metry energy. It has been proposed that the yield ratio 
of 3H/*He in HIC in wide energy range is a sensitive 
probe of nuclear symmetry energy [77-84]. However, the 
origin of 7H and *He are very complicated and the clus- 
tering process is reportedly correlated to the transport of 
isospin degree of freedom. Experimentally, in the energy 
regime of 1 GeV/u, a systematic data of the production 
of pions, protons and light clusters including 3H and He 
have been published by FOPI collaboration [17, 85], but 
in some other experiments, the 3H —? He puzzle has been 
reported [86, 87]. 

With the aim of measuring the 3H/?He ratio, the 
events generated by UrQMD model are analyzed in the 
fast simulation framework of CEEROOT. Fig.15 presents 
the 9H/*He ratio as a function of beam energy of 
208PH+28Ph system given different slope parameter of 
Esym(p), y = 0.5 (soft) and 1 (stiff), respectively. The 
charged products are filtered by the tracking detectors. 
It is shown that the yield ratio 3H/*He exhibits cer- 
tain sensitivity on Esym(p) and the difference between 
y = 0.5 and 1 is kept after the detector filtering. 
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FIG. 15. (Color online) The #H/?He yield ratio as a function 
of beam energy for reaction system 7°°Pb+?°*Pb at different 
symmetry energy parameters y = 0.5, 1.0 


In order to clarify whether the symmetry energy takes 
effect on the clustering process, one can further investi- 
gate the degree of clustering by looking at the ratio of 
the number of free protons to that of bounded in clusters 
produced in the HICs. In a general picture, the isospin 
dependent N-N cross sections, the nuclear symmetry po- 
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tential and the Coulomb interactions are convoluted ine 
forming the clusters as well as the emission of free pro-ss 
tons. Therefore, the yield of the free protons, in com-ss0 
parison to the the bounded protons, carries presumablyea1 
the information of the nuclear symmetry energy Esym(p),632 
although the evolution of the system from a compressioness 
phase at early stage to the freeze-out stage is complicatedesa 
and depends on the transport models. 635 

Fig.16 presents the excitation function of the ratio ofss. 
the free and bounded protons for the central collisionses 
208Ph+28Pb at 300 MeV/u using different y = 0.5 andes 
1.0, respectively. The yield of bounded proton counts theess 
charge of the light particles with Z < 3 in the final states 
of the reaction, including d, t, 7He and *He. The ratios 
increases gradually with the beam energy, indicating thats 
the clustering is more favored at low beam energy. More-s43 
over, it is shown that the degree of clustering, measuredesa 
by the ratio Yee / yams, can differentiate the stiffnesseas 
of Esym(p). With a softer symmetry energy (y = 0.5), 
more protons are found in the clusters. 
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FIG. 16. (Color online) The ratio of free protons to that” 


bounded in light clusters as a function of beam energy inf” 
208Pp+208Pp reactions. Two parameters of Esym(p), y = 0.5°° 
and 1.0 are compared. 


B. Radial flow and EoS of nuclear matter 


It has been realized since long time that the collective 
flow is an observable of great significance for studying 
the equation of state of nuclear matter under extreme 
conditions [88-94]. Depending on the collision geome- 
try, the collective flow has different forms. In non central 
collisions, collective motion of HIC manifests itself in di- 
rectional flow and elliptic flow, denoting the first and 
second order coefficient of the Fourier expansion of the 
azimuthal distribution of the final products, respectively 
[95]. In central collisions, the fireball produced in the col- 
lision experiences compression following by an expansion,ces 
forming the so-called radial flow [96]. It can be expectedsss 


10 


that the size of radial flow depends on the compressibil- 
ity of nuclear matter [84, 97]. Later, it is demonstrated 
by FOPI collaboration that the radial flow is sensitive to 
the incompressibility coefficient of nuclear matter [85]. 
In this subsection, the detection of radial flow with CEE 
is discussed. 

The radial flow can be extracted from the transverse 
momentum spectra of the final particles at midrapidity in 
the central collisions. Two methods have been developed 
to obtain the expansion velocity Br, one is to fit the ps 
spectra using Siemens-Rasmussen formula [98], the other 
is to extract the dependence of the mean kinetic energy 
as a function of the mass of final products based on blast 
wave model [99-101]. In both methods, the thermal 
motion characterized by the temperature T and the ex- 
pansion velocity p, are included. Here it is attempted to 
study the radial flow by the Siemens-Rasmussen formula 
[98] written as: 


1 ÆN 
pt dpdy? 


where y = 1/,/1— 62, a = 6,py/T, E and p are the 
energy and momentum of the particle in the center of 
mass reference. The radial expansion velocity 6, ther- 
mal freeze out temperature T and the constant C are the 
fitting parameters. 

Fig.17 presents the feasibility of the radial flow mea- 
surement in central collisions of ?°°Pb+?°8Pb at CEE. 
Panel (a) presents the p, spectra of proton, deuteron and 
triton. The low momentum parts are filtered by the ac- 
ceptance cut and track length cut. The curves are the 
simultaneous fitting results using formula (10). It can be 
seen that the p, spectra can be reproduced fairly well by 
Siemens-Rasmussen formula. In panel (b), the minimum 
x? distribution is plotted on the T — 6, plane. It can be 
seen that by simultaneous fitting, the two parameters can 
be determined rather well. While in individual fitting to 
the Z = 1 isotopes, the two parameters are correlated. 
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FIG. 17. (Color online) (a) The 6, and T are the fitting 
curves of the p, spectrum and the model in ?°°Pb+?°8Pb at 
0.4 GeV/u. The low transverse momentum part is eliminated 
in the fitting; (b) shows that x?/NDF is fitted as 6, and T. 


Fig.18 presents the excitation function of 6, (a) and 
T (b) reconstructed from the simulations with « = 200 
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FIG. 18. (Color online) (a) It is shown that the thermal 


freeze out temperature T increases monotonically with the 
increase of beam energy in the energy range of 0.3 GeV /u to 
0.7 GeV/u. (b) Shows radial flow 8, increases monotonically 
with the increase of beam energy in the energy range from 0.3 
GeV/u to 0.7 GeV /u. 


and 280 MeV, repsectively. The transport model sim- 
ulation reveals that different incompressibility leads to 
different freeze-out temperature T and expansion veloc- 
ity +. With a soft incompressibility constant, T is higher 
and 8r is lower, relatively, and vise versa. The difference 
is larger than 10% for most of the cases. It suggests that 
hard nuclear matter causes violent compression and ex- 
pansion, while the soft nuclear matter tends to convert 
more kinetic energy to thermal energy in the procedures 
of the collision. It is also suggested that CEE has thezos 
ability for the studies of radial flow in the collisions of, 


heavy system. e 
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C. 1 /n* yield ratio and symmetry energy 
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713 


The yield ratio of n~ /r™ produced in HIC has been 
proposed as a sensitive probe of nuclear symmetry en- 
ergy at supra-saturation densities [25, 79, 102-104]. In 
a transport picture, the pions are mainly produced by 
A resonance [33, 105]. From the isobaric model, the 
initial ratio (n~ /n+)es=(5N? + NZ)/(5Z?2 + NZ) 
(N/Z)3onsor Where N and Z are the number of neu- 
trons and protons in the reaction participating region.” 
Therefore, (m7 /7*)res is a direct measure of the isospin” 
asymmetry (N/Z)dense, Wihch is relatively higher than” 
that of the system due to the isospin fractionation ef- 
fect arising from Esym(p). In return, the yield ratio 
n/n" becomes a sensitive probe of Esym(p) [102]. In” 
a statistic model [106], the ratio t~/zt is related to 
the chemical potential of neutron and protons written7 
as nm /at ~ exp[2(un — Mp) /T], where T is the systems 
temperature, and un — Hp is the difference between the 
chemical potentials of neutrons and protons, which de-72 
pends on Esym(p) [107]. Circumstantial evidences of az 
soft Esym(p) has been found from the pion data of FOPIz 
collaboration in accordance with the isospin fractiona-730 
tion picture [20]. However, the conclusion turns to bez: 
inconclusive because of the model dependence and the7z:2 
observable dependence, since the modeling of the pionzss 
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production transport, as well as clustering, is very com- 
plicated in transport model. For this direction, one can 
refer to [36, 37, 108]. 
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FIG. 19. (Color online) The excitation function of final 


na /n* yield ratio in ?°°Pb+?°SPb reactions with different 
symmetry energy parameters y = 0.5 and 1.0, respectively. 


Therefore, pion physics is one of the focuses at CEE. 
In the current design of CEE, as shown in Fig.7(f),(g) in 
section III A, the whole phase space of pions on u® —y® 
plane is well covered, although the efficiency on azimuth 
is lost near @ = 90° and 270° in laboratory. With the cur- 
rent efficiency, the final yield of n7 /r* in ?08Pbh+°?SPb 
reactions at various beam energies are extracted and plot- 
ted in Fig.19 with y = 0.5 and 1.0, respectively. It is 
shown that the ratio m~ /r* decreases with increasing of 
the beam energy. The difference between y = 0.5 and 1 
is also more pronounced at lower beam energy, in accor- 
dance with previous predictions for Au+-Au systems [37]. 
This trend is due to the larger space-time volume of the 
fireball created in the collisions at lower energies [75]. It 
is suggested that at the favorable energy region of CEE, 
namely near 300 MeV/u , the stiffness of Esym(p) be- 
tween y = 0.5 and 1.0 can be discriminated if the sys- 
tematic error can be controlled within 10%. 


D. K? and incompressibility 


Finally, let us discuss briefly the production and detec- 
tion of kaons at CEE. Kaons are mainly produced in the 
supra-saturation density zone created in the early stage 
of the collisions, and the typical production channels are 
B+B7>B+Y+K,7+B—-Y+K. Kaon produc- 
tion is of interest in this energy region as a very clean 
probe to nuclear EOS because it experiences much less 
rescattering process in comparison to that of nucleons 
or pions. Since long time it has been pointed out that 
the yield of kaons is an effective probe for EoS of supra- 
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saturation density nuclear matter [109]. KaoS collabo-7 
ration from GSI measured the yield of Kt in C+C andro 
Au+Au systems from 800 MeV/u to 1.6 GeV/u. Withinz: 
the transport model framework, the yield ratio of Kt inm 
the two systems favors a soft nuclear matter [110, 111].7 
Compared to charged Kaons, the neutral kaon K® is notr 
affected by Coulomb interactions, so the production of 
K? carries the information of the high-density medium 
where it is produced. The yields of these neutral strange 
particles were studied in detail on the SIS, BEVALAC 
and AGS [68, 112, 113]. 

Since the HIRFL-CSR can provide proton beam of the 
maximum energy 2.8 GeV and carbon beam of the max- 
imum energy 1.2 GeV, it is also possible to conduct kaon 
experiments at CEE. Since the identification of Kt may 
suffer from the heavy background caused by protons and 
pions, it is of our interest as the first step to measure the 
neutral K? decaying to t+ and m~, which are thus easier 
to be detected. Fig.20 shows the invariant mass spec- 
trum of + and m~ in p+Ni system at « = 200 MeV. 
It is shown that, after the background are subtracted by 
mixing event method, a peak at K? mass is evidently 
seen. Totally 1 M central collision events are simulated 
and about 200 K® are reconstructed. Worth mention- 
ing, the background is not fully removed and there is an 
enhancement at low mass end, this is not due to the de- 
tector filtering, but possibly due to the residue baryon- 
baryon correlation since the majority of pions are pro- 
duced via A resonances. Compared to the Monte-Carlo 
truth tracks, the reconstruction efficiency of K° is about 
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FIG. 20. (Color online) K$ invariant mass spectrum of p+Nijo, 
system at K=200 MeV 


Fig.21 shows the yield of the reconstructed K? (mainly™ 
K?) as a function of system size with two different in- 
compressibility coefficients x = 200 and 280 MeV, re-77 
spectively. The yield of K? increases with the mass ofrs 
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the system, and exhibits significant dependence on the 
incompressibility « of the nuclear matter. As long as the 
systematic uncertainty can be controlled within 10%, the 
production of neutral Kaon K® can also be used as a 
probe of the nuclear equation of state, supporting one of 
the physical programs of CEE. 
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FIG. 21. (Color online) Variation of K$ yield in p+C, p+Ni, 
p+Ag and p+Ta reactions at beam energy of 2.8 GeV with 
two incompressibility « = 200 and 280 MeV, respectively. 


V. CONCLUSION 


The conceptual design of CEE are briefly described. 
The expected performance is studied based on the Geant 
4 simulations using UrQMD as the event generator. The 
main part of the detector system is a large-acceptance 
magnetic dipole housing the tracking detectors and the 
surrounding TOF detecotrs. TPC and MWDC are em- 
ployed to cover the midrapidity and forward rapidity re- 
gion for tracking, respectively. TOF detectors, including 
iTOF and eTOF, provide the main trigger signals con- 
structed based on advanced FPGA technology. In the 
framework of CEEROOT, built ad hoc for the researches 
and developments (R&D) of CEE, we have discussed 
the feasibility of the studies of nuclear equation of state 
through various observables, including the production of 
light clusters, 3H/*He yield ratio, radial flow, pion ra- 
tio and near-threshold Kaon production. The results of 
the simulations suggest that the aforementioned probes 
can be measured with CEE, giving the opportunity to 
the studies of properties of dense nuclear matter in the 
HIRFL-CSR energy regime. 
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